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Starting  with  the  value  of  the  temperature  gradient 
in  a  crystal  not  greater  than  10  degrees  per  centimeter,  the 
rate  of  cooling  for  NaCl  crystals  was  found  to  be  2.4  degrees 
per  minute,  for  KG1  -  2.32  degrees  per  minute,  for  KBr  - 
2.25  degrees  per  minute  and  for  KI  -  2.12  degrees  per  minu¬ 
te.  The  calculated  results  are  well  supported  by  experi¬ 
ments. 

In  the  mass  production  of  samples  from  crystals  grown 
from  a  melt  /l/,  it  was  noticed  tJ*at  in  the  cleavage  of  al¬ 
kali-halogen  crystals  in  the  direction  of  (100)  very  often  the 
splitting  surface  so  obtained  was  not  plane-parallel,  but 
wave-like.  It  can  be  assumed  that  in  the  grown  alkali-halo¬ 
gen  crystals,  as  well  as  in  metals  and  semiconductors  (3)» 
the  deformation  of  crystals  was  caused  by  thermal  stresses  due 
to  their  uneven  or  extremely  rapid  cooling. 

We  shall  consider  the  thermal  state  of  the  sodium 
chloride  crystal  grown  from  a  melt.  The  heat,  imparted  to  the 
crystal  by  fusion  is  partly  removed  by  the  cooling  water  and 
partly  dissipated  as  radiation  to  the  surrounding  medium.  We 
shall  find  the  radial  and  axial  temperature  changes  in  the 
crystal.  The  heat  transfer  equation  for  the  crystal  having 
the  form  of  a  round  cylinder  of  radius  r  in  cylindrical  coor¬ 
dinates,  is  written  in  the  form  (2) 

±L+±-*L+£L  =  0,  (I) 

dp*  p  dp  dz- 

where  $  -  is  the  flowing  radius  vector,  %  -  is  the  ordinate. 
The  limiting  conditions  for  the  solution  of  this  equation  are 
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as  follows: 

1}  When  2  =  0,  the  crystal  temperature  la  equal  to  the 
fusion  temperature. 

2)  The  amount  of  heat  supplied  to  the  JWrtal  **  «»  •*■ 
pense  of  heat  transfer  from  the  fusion  is  equal  to  the 

amount  of  heat  dissipated  by  the  crystal  surface: 


where  K  -  is  the  coefficient  of  heat  transfer,  and  A  -is  the 
coefficient  of  heat  emission.  The  °f,S?uatl°  (  } 

with  some  approximations  will  assume  the  form  (3; 


L.\ _ 
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The  ffiven  eauation  permits  one  to  determine  the  tempe- 
nature  at  aly  pSlSt  S  tMorystal.  The  analysis  of  equation 
(4)  shows  that  the  temperature  distribution  in  the  cryst, \ 
conforms  to  the  exponential  rule  in  the  direction  of  * 

and  to  the  parabolic  rule  in  the  radial  direction  of  the  a  s. 
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Fil a  i.  Axial  distribution  of  temperature  in  a 
crystal  of  NaCl  grown  from  a  melt. 

Figures  1  and  2  represent  the  results  of  calculations 
of  the  respective  axial  and  radial  distribution  of  temperature 
in  the  NaCl  crystal  with  the  assumption  that  r  -  2  centimeter, 
K  =  1  kilocalorie  /centimeter/second/ degree  (4),  and  K  - 

240  x  10-5  kilo joule/meter/second/degree  (4) *  As  can  be  seen 

f'rom  Fig.  1,  the  maximum  temperature  difference  in  th  y 
toi  in  the  radial  setting  takes  place  when  *  =  r/2  -  1 
Umeter  aud  rSLLs  M  =  124°C.  Moreover,  the  value  of  the 
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1-1  =  0.  2  —  Z  sal,  <3—2  =  2. 
4  —  z  —  3,  5  —  2  =  4,  6  —  z  — 


Pig.  2.  Radial  distribution  of  temperature  in 
a  crystal  of  NaCl  grown  from  a  melt. 


calculated  temperature  in  the  center  of  the  crystal  is  equal 
to^lS0:*,  and  at  the  surface  of  the  very  same  crystal  the 
temperature  is  equal  to  592°K» 

The  temperature  gradient  in  the  crystal  is  determi¬ 
ned  according  to  formula  (3) 

dT_  __  9  y*'  (4) 

dp  H 

where  <r  is  the  Stephan-Boltzman  constant.  The  calculation 
shows  that  when  X  =  r/2  =  1  centimeter  the  ' 

dient  is  equal  to  314  degrees  per  cent the  cm- 

f?he3gi°eat42el-.5t"SStI«e?IL?en gSflentf  the  gr4t- 
er  the  ?Sle  of .SLl  stresses  originating  In  the  crys- 
tal  under  cooling.  Since  the  temperature  is  dissimilar  in 
various  oarts  of  the  crystal,  the  mechanical  stresses  in  di- 
versHartf  of  the  crysLl  will  he  unequal.  The  central  part 
of  the  crystal  will  be  the  most  "charged  region.  The  diffe 
?ence  in  linear  dimensions  of  two  neighboring  parts  of  the 
ve?v  same  crystal  layer  in  the  direction  of  radius  r  is  pro- 
Dortional  to  the  temperature  difference  between  these  two 

that  is  Al  -  oth  T.  These  contractions  are  small, 

SH  iras-sser;  rSsI'rhSi-s  — 

ses  in  the  crystal  can  be  determined  from  the  expression 

p=ELt=EaL  T,  (5) 
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where  E  -  is  the  modulus  of  elasticity,  and  -  the  temper 
rature  coefficient.  For  NaCl,  E  =  4.2  x  10=  kilograms  per 
square  centimeter  (5),  and  oe  -  4  x  10“-'  i/per  degree  (6). 

In  the  Nad  crystal  when  &  T  =  124°C  the  mechanical 

stresses  amount  to  20  kilograms  per  square  millimeter,  that 
is,  they  exceed  its  mechanical  strength  (according  to  data 
in  /if  the  tensile  strength  for  NaCl  crystals  amounts  to 
0 .16  m  0.5  kilograms  per  square  millimeter).  Optical  stu¬ 
dies  /S/  confirm  the  presence  of  mechanical  stresses  caused 
by  the  uneven  cooling  of  the  crystal  under  growth.  Under 
the  action  of  such  significant  mechanical  stresses,  a  shift 
or  "slippage*'  of  certain  crystal  layers  relative  to  others 
in  the  crystal  facet,  parallel  to  the  surface  of  fusion,  ta¬ 
kes  place.  Defects  in  the  crystal  structure  show  up  which 
are  quite  large  in  relation  to  perfect  crystals  /9/  and  the 
density  of  which  in  the  cireotion  of  radius  r  of  the  crystal 
is  determined  from  the  relationship  /3 / : 

w— n.JdL.  <6) 

a0  Ap 

where  AT/^f  is  the  radial  temperature  difference.  Recent 
optical  and  electrical  measurements  /8,  9/  confirmed  the  hy¬ 
pothesis  regarding  the  existence  of  structural  micrcfractures 
in  deformed  crystals.  The  lattice  defects  change  the  elec¬ 
tron  mobility  thus  contributing  to  the  accumulation  of  a 
space  charge  which  in  turn  influences  the  crystallizability 
of  matter.  For  the  determination  of  the  density  of  these 
defects,  we  shall  assume  that  A  T /&f  =  10  degrees  per  centi¬ 
meter  /10,  11/.  X-ray  studies  /12/  with  NaCl  have  shown  that 
there  are  no  mechanical  stresses  operating  when  f  =  10 

degrees  per  centimeter.  The  calculation  for  our  case  gives 
N  5  1.75  x  10 ^  ^/square  centimeter.  Xhe  rise  of  radial  tem¬ 
perature  differences  results  in  an  Increase  In  the  number  of 
lattice  infractions.  Therefore,  in  the  growing  of  crystals 
from  fusion,  a  significant  role  is  played  by  the  rate  of 
crystal  cooling  during  the  growth  process,  which  can  be  de¬ 
termined  according  to  formula  /3/ 

(7) 

dt  3  cr  A  p 

where  c  -  is  the  specific  heat  of  the  crystal.  For  NaCl, 
c  =  0.24  calories/gram  degree  /13/.  The  cooling  of  the  crys¬ 
tal  should  be  carried  out  slowly  so  that  the  radial  tempera¬ 
ture  gradient  be  kept  at  a  minimum.  Expressing  it  differen¬ 
tly,  the  isothermal  surfaces  in  the  crystal  should  be  paral¬ 
lel  to  the  surface  of  fusion.  The  cooling  rate  of  the  NaCl 


4 


crystal,  as  determined  by  the  use  of  formula  /6/  with  ^T/- 
A?  =  10  decrees  per  centimeter  amounts  to  dT/at  =  2 J16” 
g/ees  per  minute,  that  is,  of  the  same  orger  oj.  magnitude  as 

in  /12,l4/. 

The  crystal  will  grow  only  In  case  the  rate  of  growth 

will  be  higher  than  its  rate  of  cooling*  It  Is  known  - 

the  velocity  *Vh  ,  with  which  pure  crystals  grow  from  fusion, 
can  be  determined  approximately  according  to  formula  /15/ 


-T”'  Wl 

R»Tn  Ac 

where  D  is  the  diffusion  coefficient,  which  for  NaCl  is  equal 
to  0.14  square  centimeters  per  second;  S,tne  molar  entropy 
of  fusion  is  equal  to  38  calories  per  mole  degree,  Az  is  the 
distance  in  the  direction  of  the  crystal  s  expansion,  in 
which  the  crystallization  takes  place  (assumed  to  be  equal  to 
5  x  10*2  centimeter  /l5/,  R0  is  the  universal  gas  constant, 
and  ATr  -  is  the  drop  in  temperature  at  the  crystallization 
boundary  (at  the  distance  Az  )  •  From  Fig.  1,  for  Ax  -  5  x 
10-2  centimeter  with  AT/hy  =  10  degrees  per  centimeter,  we 
find  Aik  =5  degrees,  and  the  rate  of  growth  "V*  - 
centimeters  per  second,  that  is,  a  significantly  higher  cool 

ing  rate. 

For  other  alkali- halogen  crystals  an  analogous  calcula¬ 
tion  with  several  additional  approximations  can  be  carried  out. 
We  were  unsuccessful  in  an  attempt  to  determine  the  heat  emis¬ 
sion  coefficient  K  of  other  crystals.  For  NaGl  the  value  oi 
the  ratio  X  /K  is  0.175.  l/centimeter.  This  ratio,  roughly 
speaking;,  will  be  constant,  since  the  very  same  coefficients 
and  K  depend  to  a  small  extent  upon  temperature,  and  have 
values  of  the  same  order  of  magnitude.  It  can  be  assumed 
that  the  ratio  of  these  coefficients  is  axso  equal  to  o.lfb. 
l/centimeter  for  the  rest  of  the  alkali- halogen  crystals. 

We  consider  that  the  error  introduced  with  this  as¬ 
sumption  is  small  and  is  found  to  be  within  the  limits  of 
precision  of  the  performed  calculations.  The  values  for  th 
rate  of  cooling  as  obtained  from  calculation  |re  as  follows^ 
for  KOI  -  2.32,  for  KBr  -  2.25,  and  for  KI  -  2.12  degrees  per 

minute. 

The  deduced  calculation  shows  that  the  extremely  ra¬ 
pid  cooling  of  alkali-halogen  crystals,  grown  from  a.  melt,  can 
be  the  cause  of  significant  mechanical  stresses  the  occur¬ 

rence  of  structural  defects  in  crystals.  The  values  found 
for*  crystal  growth  rate  are  supported  by  data  of  other  authors 


/12.  14/,  and  as  the  performed  experiments  demonstrate, 

'they  can  be  utilized  In  the  practice  of  growing  aljcaxi-halo- 
gen  crystals  from  a  melt.  The  rate  of  crystal  cooling  ob¬ 
tained  in  the  given  calculations  can  he  assumed  to  ® 

maximum  cooling  and  heating  rates  of  samples  from  alkali-ft*- 

iflstAjaauia  .  „  ,  .  +  'rtenftririfin- 


Bibliography 

1,  Kvropoulos,  3*  Zelt.  f.  anorg.  u.  alls*  Ghera,  154,  308, 

1926 

2,  Lykov,  A.  V*,  Teoriya  teploprovodnosti  (Theory  of  Heat 

Transfer) .  GrIITTL,  1952* 

3,  '  Billig,  E.  Proc.  Roy,  3oc»,  235,  37,  1956, 

4„  Val'ter,  A,  F.,  Inge,  L.  D*  ZhTP  (Journal  of  Technical 
Physics),  8,  295,  1958* 

5.  Dorfman,  Ya.  <J.,  Frish,  S.  E.  Editors?  Sh.  fizicheskikh 
Constant  (Compendium  of  Physical  Constants),  OHTx, 

1937. 

6  Zluze,  V.  P.  PAN  SSSR  (Reports  of  the  Academy  of  Scien¬ 
ces  USSR),  99,  711,  1954. 

7 „  Kuznetsov,  V*  D.  Fizika  t'verdogo  tela  (Physics  of  So¬ 
lids}*  Izd.  “Krasnoye  anamya  ,  Tomsk,  1941, 

8,  Billig,  E.  Brit,  J*  Appl.  Phys.,  7,  375,  1956. 

9*  Cottrell,  A*  H*  Nature,  178,  1090,  1955. 

10.  Strong,  J.  Phys.  Rev.,  36,  1.663,  1930. 

11.  Sears,  G*  M. ,  Coleman,  R.  Ve  J.  Chea.  Phys,,  25,  685, 

i  a  r-;A 


12*  Jones,  D*  A. ,  Smith,  T.  Proc.  Phys,  Soc.,  B  69,  8?8,  1956. 

13.  Baehinskiy,  A.  I*,  Putilov,  V*  V,,  Suvorov,  N*  P*  bpra- 
vochnlk  po  flzike  (Handbook  of  Physics).  Uchpedgiz, 
1951. 

14*  Vermeee,  J®  Physic a,  22,  1257,  1956* 


6 


15*  Hill ig,  W,  B.,  and  Turnbull,  De  J.  Ghem,  Pbys.»  24,; 
914,  1956* 


Submitted  for  pu¬ 
blication,  6  Novem¬ 
ber  1957* 


.END 


1471 


Tomsk  Polytechnic  Insti¬ 
tute  imenl  !&*  kirov* 


n 


FOR  REASONS  05*’  SPEED  AND  ECONOMY 
THIS  REPORT  HAS  BEEN  REPRODUCED 
ELECTRONICALLY  DIRECTLY  FROM  OUR 
CONTRACTOR'S  TYPESCRIPT 


This  publication  was  prepared  under  contract  to  the 
UNITED  STATES  JOINT  PUBLICATIONS  RESEARCH  SERVICE 
a  federal  government  organization  established 
to  service  the  translation  and  research  needs 
of  the  various  government  departments. 


